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When subjected to a horizontal temperature difference, a fluid layer with a 
free surface becomes unstable and hydrothermal waves develop in the bulk. Such 
a system is modelized by two coupled amplitude equations of the one- dimensional, 
r^' ' complex, cubic Ginzburg-Landau type. By transposing the method developed for 

(*— ■) , one CGL3 equation, we obtain several new exact solutions expressed by closed form, 

f"*) ' singlevalued, analytic expressions. Some of them are the analogue of the famous 

amplitude hole solution of Bekki and Nozaki. 



1. Introduction 

Consider a thin fluid layer with a free surface, subjected to a horizontal difference of 
temperature between its vertical edges. When the temperature difference T%— Tj is small, 



> 

•i-H 

X. 

the basic flow is just one long, narrow convection cell. If one increases T% — T\, this basic 
flow becomes unstable via a supercritical Hopf bifurcation, and traveling hydrothermal 
waves appear (I], ||. Two waves develop (see Fig. f in Ref. Q), and their complex 
amplitudes (A + , A_ ) evolve in time according to 0] 

d t A± = rA ± T vd x A ± + (f + ia)AA ± - (f + i/3)(|^l ± | 2 + 7^tI 2 )^± ; (1) 

in which r is real and proportional to T% — Ti, while (w,a,/3,7) are real parameters and 
A is the Laplacian. 

These are two coupled nonlinear partial differential equations of a very generic type, 
namely the cubic complex Ginzburg-Landau one (CGL3). 

TO APPEAR IN REPORTS ON MATHEMATICAL PHYSICS. S99/054 

[1] 
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2. The problem for two coupled one-dimensional CGL3 equations 

Let us for simplification consider one space dimension. 

The scalar one-dimensional CGL3 equation for A(x, t) is defined as 

E = iA t + P A xx + q\A\ 2 A-i 1 A = 0, pq^O, {A,p,q)(=C, 7 eK, (2) 

with p, q, 7 constant. The review Gj recalls the physical phenomena (pattern formation, 
superconductivity, nonlinear optics, . . . ) modelized by this equation and summarizes the 
known exact results. 

Two coupled one-dimensional CGL3 equations are similarly defined as 

f iA t + ivA x + P A XX + q(\A\ 2 + 6\B\ 2 )A - ijA = 0, , , 

" \ iBt - ivB x + P B XX + q(\B\ 2 + S\A\ 2 )B - ijB = 0, l ' 

in which the coupling parameter 8 is a complex constant and v is the group velocity. These 
describe for instance the amplitudes of two lasers [pj , or the spatiotemporal intermittency 
[j, 0] , or hydrothermal waves g] . 

The problem which we address is to find exact solutions, i. e. analytic expressions 
written in closed form, able to describe the observed topological structures : pulses, 
fronts, holes, . . . The closed form which we require excludes any perturbation process, 
which would involve infinite series generically unable to globally represent the solution. 

More specifically, we restrict to the search for the most general solitary wave 

A(x,t)=a(t)e< UJAt + ^ A ^\ B(x,t) = b(0e< ujBt + ^ B ^\ f = x - ct, (4) 

in which (o, 6, if a, <^b) are functions of the reduced independent variable £, and we only 
consider the strict CGL3 case Im(p/q) ^ 0. 

3. Analyticity restricts the number of free constants 

The chaotic nature of CGL3 restricts the number of arbitrary constants the solitary 
wave (0) may depend on. Let us explain how on the Kuramoto-Sivashinsky equation, for 
which the solitary wave reduction is the ordinary differential equation (ODE) 

vu" + iiu + u" /2 + K = 0, (/i, v) real constants, (5) 

in which K is an integration constant. Denoting \ = x — xq, the general solution is 
locally represented as ||, || 

u(a:o, c+,c_) = 12(V?T 3 {Taylor(x) 

+ c + x( 13 + J ^/ 2 (. • •) + C -^ 13 ~ ^/ 2 • ■)}, (6) 

in which (xq, c+,cJ) are arbitrary constants, Taylor a converging Taylor series of %, and 
the two (. . .) highly multivalued series depending on \ with only positive integral powers 
of (c + , C_) . The presence of chaos is linked to the fact that the series (a) is what Painleve 
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calls "une solution illusoire" , and the only way to possibly recover analyticity is to require 
c + = c_ =0, restricting the illusory general solution to what can be called the general 
analytic solution, depending on 3 — 2 = 1 arbitrary constant, namely xq. 

The above two irrational exponents are computed as follows fllpf , in the example of 
the ordinary differential equation (|^). 

In a first step, one looks for a singular dominant behaviour u ~ uqx p (with uq nonzero 
and p not a positive integer) when x — > xq. The term u' is less singular than u'" , the 
term K is regular and cannot contribute, and so the dominant behaviour is governed by 
vvl" + u 2 /2, which contributes the terms 

MP - 1)(P - 2 WX^ 3 + (1/2) (u oX p ) 2 , (7) 

and generates the two conditions 

p - 3 = 2p, MP - l)(p - 2)«o + (l/2)u§ = 0, «o ?* 0. (8) 

Their unique solution is p — —3, u a = 120^, and the common value of the two powers is 
q=p — 3 = 2p = —6. 

In a second step, one builds the linearized equation near the solution which behaves 
like uqX p 

(i/flg + u aX p )w = 0. (9) 

This linear ODE for w has a Fuchsian singularity near \ = (see the textbook |pT|), and 
one then computes its Fuchs indices, i. e. the roots j of the algebraic equation 

\im X - 3 - q (vd 3 x +u x p )x 3+P (10) 

= v{j - 3)(j - 4)(j - 5) + 120^ = u(j + l)(i 2 - 13j + 60) (11) 

= "(J + 1) [J S \\J o = °- ( 12 ) 




A necessary condition for analyticity is that the two irrational Fuchs indices do not 
contribute, and this restricts to 3 — 2 = 1 the number of arbitrary constants which the 
general analytic solution may depend on. 

Let us now compute the similar quantities uo,p and the set of Fuchs indices for the 
(much more difficult) coupled CGL3 system. 

4. Leading order of two coupled CGL3 

One easily checks that |A| and \B\ generically behave like simple poles. Let us then 
denote the dominant behaviour of the four fields (A, A, B, B) as [jjjj 

A - a oX - 1+ia , A ~ a oX ^ ta , B - 6 x^ 1+ ^, B - &oX~ 1_i/3 , (13) 

in which (ao,6o) are complex constants, (a,/3) real constants, all of them to be deter- 
mined. The analogue of p — 3 = 2p is now just the identity, and the analogue of the 
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second equation in (|8j) is the set of four equations (the time-derivative term and the 7 
term are not dominant) 

p(-l + ia){-2 + ia) + q(a 2 + 5b 2 ) = 0, (14) 

p(-l - ia)(-2 - ia) + q{a 2 + ~5b 2 ) = 0, (15) 

p(-l + t/9)(-2 + z/3) + #2 + <5a 2 ) = 0, (16) 

p(-l - i(3)(-2 - if3) + q(b 2 + 6a 2 ) = 0. (17) 

with a 2 — \ao\ 2 ,b 2 = \bo\ 2 . These are four nonlinear algebraic equations. The methodo- 
logical error not to be done would be to solve them for the four real unknowns (a 2 , b 2 ,a, (3) 
as real expressions in the three complex parameters (p,q,S), for this would generate 
algebraic (i. e. not rational) expressions and mess up all subsequent computations. On 
the contrary, one should make no distinction between the unknowns and the parameters, 
and consider this apparently nonlinear system as a linear system on the complex field 
C for a carefully chosen subset of four unknowns among the ten equivalent variables 
(a 2 , b 2 , ia, i/3,p, q,p, q, S, 5) so as to avoid the introduction of any algebraic quantity. Such 



a subset is (a 2 , b 2 ,p, q). First, the subsystem (14), (|16|), linear inhomogeneous in (02, b 2 ), 
is solved on C as 

a 2 = q{1 P _ s2) (-(! - *«)(2 - ia) + 6(1 - ij3)(2 - i0)) , (18) 

b 2 = q{1 P _ 52) (-C 1 - W - i0) + 5(1 - ia)(2 - ia)) , (19) 

if one excludes the nongeneric case S 2 = 1, which is left to the reader. Second, the 
subsystem (|15|), (0) is linear homogeneous in (p,q), therefore its Jacobian must vanish 
(from now on, one restricts to S = 5) 

p B{p,q) 

and the subsystem is solved on C as 

p = Kp(\ - ia)(2 - ia), q = Kq(l + ia)(2 + ia), (21) 

in which K is an irrelevant arbitrary nonzero complex constant. 

The constraint (EQ) between the real variables (a, (3), which describe the argument of 
(A, B), defines two mutually exclusive cases 

1. a ~ (3, a priori not so interesting since it implies a 2 — b 2 and contains the reduction 
A = B to one CGL3 equation, 

2. a(3 = — 2, which forbids a — (3 and describes a truly coupled situation. 
For reference, the resolution on 1Z provides the algebraic solution 

a >= w^m [{1 - 6)dr+{1 - eS)A] ' (22) 
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2|#(l-tf»)d? 

2^ + A ^ 2,/,. 



[(I - 5)d r + (e - S)A], (23) 



r^-(d r +A), p=JL(d r + eA), (24) 



P(j) 



- = d r - id h A 2 = d 2 r + (8/9)d 2 , e 2 = 1, (25) 

which defines four families of movable singularities pC] for the dominant behaviour. 

5. Fuchs indices 

The indicial equation is now the determinant (lfj of the fourth order matrix 

/ (2a W> + Sb b )q ga 2 , _ Sqa ba Sqa b \ 

qa^ 2 (2a a^ + 5b b a )q _Sqa^b Q 6qa b 

5qaob Sqa Q b (2b b + 5a ao)q qb 2 , 

V Sqa^bo Sqa a b qb {2b b + 5a W>)q J 

+ diag(p(j - 1 + ia){j - 2 + ia),p(j - 1 - ia)(j - 2 - ia), 

p(j - 1 + i/3)(j - 2 + ip),p(j - 1 - i/3)(j - 2 - t/3)). (26) 

In the limit a = f3 = 0, it admits the symmetry Vj P(j) = P(3 — j). After substitution of 
the values of (a2,b 2 ,p,q) as given by (|l8|), (|19|), (|2l|), the variables (p,q,K) are factored 
out and the determinant is a polynomial of degree eight in j with coefficients polynomial 
in the real variables (a, (3,5). Three of its zeros j are already known : — 1 as always in 
this kind of indicial equation, twice, which correspond to the phase invariance of A and 
B. Our interest here is to count the number of irrational indices, as done in Section |3.| 
Let us split the discussion according to the two subcases (|20|). 
When a = (3, the indicial equation factorizes as 

detP(j) = (j + l)(j 2 -7j+6a 2 + 12)j 2 x (27) 

[U + 1)(J 2 - 7j + 6a 2 + 12)(1 + S) + 4*((2 - a 2 )j - 3(2 + a 2 ))] .(28) 

If the coupling 5 vanishes, one recovers as expected the square of the indicial equation of 
one CGL3 |jl|. For generic values of the fixed parameters (p,q,5), all indices, excepted 
(—1,0,0), are irrational. 

When a[3 = —2, the remaining fifth degree polynomial does not factorize, and the 
indicial equation is 

det P(j) = (j + l)j 2 [(1 - S 2 )(j 5 - 13j 4 ) + 40 terms] = 0. (29) 

Again, there generically exist five irrational indices. Some particular values are 

a = -/3, det = (j + l)j 2 (j 2 - 7j + 24)((1 - <5)(j 3 - 6j 2 + 17 j) + 24(1 + 5)), (30) 

6 = 1, det = {j + l)j 2 (a 2 j 2 -(<i + 3a 2 + a 4 )j + 3{2 + a 2 ) 2 ), (31) 

d = -l, det = (j + l)j 2 (a 2 (/-7j 2 ) + 3(2 + a 2 ) 2 (j-3)). (32) 
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6. Counting of the free constants 

Let us count the precise number of arbitrary constants the general analytic solution 
for the solitary wave reduction (|J) may depend on. The reduction introduces the three 
arbitrary constants (c,w^,Wb) and does not change the differential order eight of the 
PDE system. From the eleven possible arbitrary constants, one must subtract 

1. the irrelevant origin £o of £ (Fuchs index —1), which represents the invariance under 
a space translation, 

2. the irrelevant origin ipo of each phase (Fuchs indices (0,0)), which represents the 
invariance under a phase shift, 

3. and the number of irrational Fuchs indices, generically five. 
The result of this counting is summarized in Table [l| 

Table 1: Maximal number of arbitrary constants in the solitary wave of one and two 
CGL3 equations. 



CGL3 


Reduction var 


Diff. order 


Irrelevant 


Irrational Fuchs indices 


Max arbitrary 


1 


C, U! 


4 


£o,¥>o 


2 


2 


2 


C,0JA,U>B 


8 


£o, </U,o, Vb.o 


5 


3 



Another counting, for the various possible topological structures, is made in Ref. p4| . 

7. Hints from the integrable limit 

In order to have an idea of the class of expressions among which to search for exact 
solutions, let us consider the integrable limit of our system, i. e. the unique case p real, 
q real, 7 = and, for two coupled equations, v = 0, S = ±1 |lq, |lq ]. 

For one equation, this is the famous nonlinear Schrodinger equation. The indices 
(— 1, 0, 3, 4) are all integral, and the general solution for the solitary wave 



\Af - 


= a 2 (p(£) - e ), 




(33) 


A 


= ^+^-£ + (l/2)^-(4e3- 
= p r ((c/(2p r )) 2 + 3e ), 


ac a) f ^ 


(34) 
(35) 


arg/i 

bJ = 


5iU ga) J p(0-eo' 



in which p is the Weierstrass elliptic function, depends on the four arbitrary constants 
(c, oj, 32,33)- 

For two equations, the eight Fuchs indices (—1,0,0,0,3,3,3,4) are integral and the 
general solution for the solitary wave, which depends on eight arbitrary constants, is not 
yet known. Recently Porubov and Parker |17| found a six-parameter elliptic solution 



\A\ 2 =a 2 (p(0- ei ), 
\B\ 2 =b 2 (p(t;)-e 2 ), 



(36) 
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in which (02,62) are those of Section [h|, ei,e2 are constant, and also, because S 2 = 1, a 
five-parameter elliptic solution 



I^I 2 ^i(p 2 (£) + 5ip(£) + ^i), 



2 _ .1 , .2,, , , „ ,<-,,/-., (37) 



in which A&, Bk, C& are constant. 

When the elliptic function p degenerates to a trigonometric function, the six-parame- 



ter solution ( 36 ) includes all the various front-type or pulse- type solitary waves [R8L [L6| 



8. What was known for one and two CGL3 

For one CGL3, one knows only four particular solutions of the reduction £ = x — ct 
with a zero codimension, i. e. without constraint on the fixed parameters (p, q, 7). In all 
of them \A\ 2 is a degree- two polynomial in tanhfc£, i. e. a trigonometric degeneracy of 
the integrable case (|3^). These are qualitatively (K or Kj denotes a real constant, and 
for brevity we omit the argument k£ of the functions) 

1. a pulse or solitary wave (lJJ 

\A\ 2 = a 2 sech 2 , c=0, (38) 

2. a front or shock [£0| 

\A\ 2 = a 2 (tanh±l) 2 , (39) 

3. a source or propagating hole J21| 

|A| 2 = a 2 ((tanh+#i) 2 +i£f) , c arbitrary, (40) 

4. an unbounded solution [M 

|A| 2 = aa^an 2 -!-^ 2 ), c = 0. (41) 

Only the propagating hole depends on one arbitrary constant (the velocity c). The 
missing solution should be, as indicated in Table [y and explained in detail in Ref. [J22| , 
an extrapolation of this propagating hole to one more arbitrary constant. 

Several other solutions have been "observed" in numerical simulations or in analytic 
perturbations and quite certainly correspond to unknown analytic solutions, these are 

1. a propagating pulse |23], i. e. some extrapolation of ( J38| ) to an arbitrary velocity, 

2. a homoclinic propagating hole [E4J, i. e. a solution with the shape of ( f40| ) but with 
the same limit of \A\ at both infinities £ — > ±oo, 

3. various other topological structures p5[ . 

For two coupled CGL3 equations, the only analytic attempt we are aware of ]26| 
results into equal phases for A and B, i. e. essentially into the uncoupled situation. The 
result to be found should, according to Table m, depend on three arbitrary constants. 
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9. Search for singlevalued global solutions 

The difficulty comes from the fact that the natural physical variables (A, A) or 
(ReA,ImA) (it is sufficient for the moment to consider one CGL3) are always locally 
multivalued |13| , hence not adapted to the search for singlevalued exact solutions. A 
more detailed study |22f uncovers the best representation for this purpose, namely a 
complex modulus Z and a real argument uniquely defined by 

A = Ze ie , A = Ze~ m . (42) 

For each of the families enumerated in Section [l], if one excludes the contribution of the 
irrational Fuchs indices, the three fields {Z, Z,grad0) are locally singlevalued and they 
behave like simple poles. 

The physical variables {\A\ 2 , gradarg^4) also have this nice property of being locally 
singlevalued (they respectively behave like a double pole and a simple pole), but they 
are not as elementary as (Z, Z, grad0). 

In order to obtain singlevalued exact solutions, one looks for representations, if they 
exist, of these locally singlevalued fields by Laurent series which terminate, thus ensuring 
ipso facto their closed form. This is the famous "truncation method" initiated by Weiss, 
Tabor, and Carnevale [M, the latest version of which is detailed in Ref. p8| . 

For two coupled CGL3 equations, this method can be applied either to the two couples 
(|A| 2 ,gradarg A) or to the two couples (Z,Z,grad0). The first set can in principle, 
despite the much more involved computations, yield more solutions than the second set, 
but it should be noted that, for one CGL3, all known solutions are found with the second 
set. In order to shorten computations, let us handle the two couples (Z, Z,gradO). 

10. How to minimize the computations 

Quite similarly to Section J4], the natural methodology (to solve for real variables 
as real expressions) is not the one to follow, let us explain it on the truncation which 
provides the propagating hole solution ([h]) of one CGL3, since the two-CGL3 case will 
follow the same method. 

The assumption is (2^] 

f Z = a Q ( X ~ 1 +X + iY), 

Z = a ( X - 1 +X-iY), 

Q = ojt + aLog^ + K^, { ' 

{ (Log^)' = X-\x' = l-(fc74)x 2 , 

in which \ an d ^ are functions of £ = x — ct, (ao,a) are the constants of Section fO 
(uj,X,Y,K,k 2 ) are real constants. After elimination of any derivative of Log?/' and x, 
the lhs E of the CGL3 becomes a Laurent series which terminates 

3 

i=o 
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and one has to solve the four complex (eight real) equations Ej = in the eight real 
unknowns (02, a, u), X, Y, K, c, k 2 ), the two complex parameters {p,q), and the real pa- 
rameter 7. If there exists a solution, the elementary building block functions evaluate 
to 



X 



k , *£ , 
-tanhy, -0 



cosh 



H 



(45) 



The good methodology is again to select, among the eleven complex variables consid- 
ered as equivalent, four variables which make the system a linear one of Cramer type. In 
this example, one proceeds as follows [p2l. The system (Eq, Ei,E<i) is of Cramer type in 
(a2,K,uj), and after its resolution the last equation £3 is independent of (p,q,j,c) and 
factorizes into a product of linear factors 



E, 



[k 2 - 4(X + iY) 2 ](aY - 2X) = 0. 



(46) 



Let us apply this to the coupled CGL3 system. 



11. Double holes 

The analogue of the hole solution of one CGL3 [|l] is searched with the assumption 

A = aoix- 1 +X a + iY a )e + ^ UJ - t + « L °gV- + K a £,) ^ 
A = Mx~ l +X a - ^y a )e- ^ (^ i + "Log V + K a £) ^ 
B = b^- 1 +X b + iY b )e+ l ^bt + /3Log^ + Kg) ^ 
[ B = boix' 1 +X b - iY b )e~ l ^ ht + ^ L °S^ + Kb Z\ 



(47) 



with the same definition for (Logip)' and x' as m previous Section. These are eight com- 
plex equations in fourteen real unknowns (02, 62, ck, (3, u> a , u> b , X a , Y a , X b , Y b , K a , K b , c, k 2 ), 
two complex parameters (p,q), and three real parameters (v, 7, S). After solving the six 
equations (Eq,Ei,E2) as a Cramer system on C in (02, 62, K a , K b ,u> a , u>b), the last two 
equations E3 factorize into a product of linear equations: 



[k 2 - A(X a + iY a ) 2 ] F(Y a , X a , X b , a, 0) = 0, 
[k 2 - A{X b + iY b ) 2 ] F(Y h , X b , X a , [3, a) = 0, 
F(Y, X a , X b , a, 0) = aY- (X a + X b ) 

v ^-ia-6(l-iP) + 6 2 (l~ia) 2 
(1 — d z )(l — 101) 

The bifurcations successively encountered during the resolution on 1Z are 



3=0 

.7 = 1 
3 = 2 
3=3 



(j3 = a) or (a/3 = -2), 
Pi zero or nonzero, 
(a — j3)p r zero or nonzero, 
the three subcases of (48)-(|l). 



(48) 
(49) 

(50) 



(51) 
(52) 
(53) 
(54) 
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The solutions p r Pi 7^ which do not reduce to a mere product of two CGL are 

X b = X a , aY a = f3Y b = 2X a , v = 0, a/3 = -2, c = & P r l \p\ 2 X a , 

k 2 = -4(3 + (a + 0)pi/pr)(3 + {a + /3) 2 5/{l + S))X 2 , 
X b = X a , Y a = Y b = 0, v = 0, k 2 = 4X 2 , a/3 = -2, c = &pT l \p\ 2 X a , 
X b = -X a , Y a = Y b = 0, c = 0, k 2 = 4X 2 a , a = {3 
X b = -X a , Y a = Y b = 0, k 2 = AX 2 , a/3 = -2, 



-2{a 2 -f3 2 )pi l \p\ 2 X a 5/(l-5 2 ), 



v = -2((1 + g)(3 + S) + (a + P) 2 S) P r 1 \p\ 2 X a /(l - S 2 ), 
X b = X a , Y a = 0,(3Y b = 2X a , v = 0, c= Qp^pfXa, k 2 = 4X 2 , 
Y a = Q,/3Y b = 2X b , k 2 = 4X 2 , a(3 = -2, 

c=(X a + 5X b )\p\ 2 / Pt ,v = (X b - X a )\p\ 2 / Pl , 



(55) 
(56) 
(57) 



(58) 
(59) 



(60) 



with the respective codimensions 1,1, 1,2, 2, 2. The fixed constraints are respectively 



v = 0, 
v = 0, 



|p| 2 (3-<5) 2 7- Pi(S,p r ,pi,a)piV 2 = 0, 

|p| 2 ((l + <5)(3 + 6) + {a + /3) 2 <5) 2 7 - P 2 (S,p r ,Pi, a + (3) Pl v 2 = 0, 

P 3 (5,p r ,pi,a + /?) = 0, 
v = 0, P4,(6,p r ,pi,a,/3) = 0, 



(61) 

(62) 
(63) 

(64) 
(65) 



a" 



+ 25 = 0, P 5 (p r ,pi,a)j 2 + P 6 (p r ,pi,a)jv 2 + P 7 (p r ,pi,a)v 4 = 0, (66) 



in which P n denotes a polynomial of its arguments. 

These "double hole" solutions depend on no arbitrary at all (as compared to the max- 
imal value of three, see Section pj), they are all heteroclinic and can never be homoclinic. 

12. Double pulse solution 

A double pulse similar to ( p8[) is searched with the assumption 

A = {aoOT 1 +X a + zF Q ) + {k 2 /A)a lX } e+*( w »* + ^ + ^0, 
A = {MX- 1 +X a - iY a ) + (k 2 /4)a lX ) e~^ f + ^ + *a0 , 
B = {bo{x- 1 +X b + iY b ) + (fc 2 /4)6 lX } e +i ^t + Vb + ^0, (67) 

B = {MX- 1 +X b - iY b ) + (fc 2 /4)6 lX } e ~ J M + Vb + ffbO, 
[ ^ = ax" 1 + (fc 2 /4)«ix, y' B - ftT 1 + (fc 2 /4)/3iX, 

in which ao and a\ are two different roots of the leading order equation, so as to produce 
a sech according to the elementary identities J22j1 



tanhz 



1 



tanhz 



-2isech 



2z + t- 



tanh z 



1 



tanhz 



2 tanh 



2* + ,- 



(68) 
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This generates fourteen complex equations in the eighteen real unknowns 
(a ,b Q ,ai,bi,a,/3, ai, f3i,u) a ,u> b ,X a ,Y a , X b , Y b ,K ai K b ,c, k 2 ), the two complex parame- 
ters (p, q), and the three real parameters (u,7, 5). Among the eight cases of signs in (p4|), 
let us restrict to the two which allow a 2 , = a 2 and 6q = b\ ; their sech-type solution is 

«i = -ao, Kl 2 = a 2 , ai = a, X a = 0, Y a = 0, K a =(c- v)p r / (2\p\ 2 ) , (69) 

h = -b Q , \b Q \ 2 = b 2 , A = (3, X b = 0, y b = 0, K b = (c + V K/(2|p| 2 ), (70) 

Cft = 0, «p( = 0, (a - /3)[(a + (3) Pl + 2 Pr ] = 0, (71) 

7 = ((1 - (a 2 + {3 2 )/2) Pt ~(a + f3)p r )k 2 . (72) 

The two moduli are proportional to sech kx but the two phases are different if a =/= /3. 
The codimension is either one (case a = /3) or two (case a(3 — — 2), and, if p; = 0, the 
two velocities are arbitrary. 

13. Conclusion 

From the singularity structure, we have found several exact solutions to the coupled 
CGL3 system. At the time being, these are the only known analytic solutions to this 
system. Their main features are 

1. their physical relevance, since they represent structures observed in experiments 

§01. 

2. their hole or bell-shaped profile characterizing their nonlinearity (these are not 
plane waves), 

3. their true coupling. 

However, a drawback is their nonzero codimension. 

In particular, Bekki-Nozaki holes have been experimentally observed in a rectangular 
geometry with a nonzero group velocity v |2J]. 

The gap between the found results and the expected results described in Section J6] 
is much more important for the two-CGL3 system (three missing arbitrary constants, 
nonzero codimension) than for one CGL3 (one missing arbitrary constant, zero codimen- 
sion). Therefore the effort should be put on filling the one-CGL3 gap, before tackling 
again the two-CGL3 system. 
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